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Abstract 

o 

The ability of Gaussian noise to induce ordered states in dynamical systems is here presented in an overview of the 
main stochastic mechanisms able to generate spatial patterns. These mechanisms involve: (i) a deterministic local 
dynamics term, accounting for the local rate of variation of the field variable, (ii) a noise component (additive or mul- 
tiplicative) accounting for the unavoidable environmental disturbances, and (iii) a linear spatial coupling component, 
which provides spatial coherence and takes into account diffusion mechanisms. We investigate these dynamics using 
analytical tools, such as mean-field theory, linear stability analysis and structure function analysis, and use numerical 
simulations to confirm these analytical results. 
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1. Introduction 

Spatial patterns are widely present in different natural dynamical systems. Their occurrence has been studied for 
quite a long time with applications to different fields, including for example hydrodynamic systems (e.g. Rayleigh- 
Benard convection [1, 2|) and biochemical and neural systems (see, for instance, ||3] |4)). In particular, a number 
of environmental processes are known for their ability to develop highly organized spatial features. For example, 
remarkable degrees of coherence can be found in the spatial distribution of dryland and riparian vegetation |5 6,[7]|8)), 
river channels J9] [TD] [TT] , coastlines ifTZlfHl . sand ripples and dunes 1141 . Figure[T]shows an example of natural spatial 
patterns around the world. These patterns exhibit amazing regular configurations. Found over areas of up to several 
square kilometers, they can occur on different soils and with a broad variety of vegetation species and life forms 
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The study of patterns can offer useful information on the underlying processes causing possible changes in the 
system. In recent years, several authors have investigated the mechanisms of pattern formation in nature, and their 
response to changes in environmental conditions or disturbance regime. For example, in the case of landscape ecology, 
these studies have related vegetation patterns to the underlying eco-hydrological processes EES US), the nature of 

T^-j- the interactions among plant individuals ll6l lT9ll , and the landscape's susceptibility to desertification under different 

("\| climate drivers and management conditions l2Tll22ll . 

ly-^ Deterministic mechanisms of pattern formation have been widely studied ll23l l24l l25l with a number of applica- 

tions to environmental processes l26l l8l l6l l27l l28l l29l [30l [3D . Stochastic models have only been developed more 
recently Il32l[33ll . They explain pattern formation as a noise-induced effect in the sense that patterns can emerge as a 
consequence of the randomness of the system's fluctuations. These random drivers have often been related Il34l[35l 
to a symmetry-breaking instability. They destabilize a homogeneous (and, thus, symmetric) state of the system and 
determine a transition to an ordered phase, which exhibits a degree of spatial organization. In the thermodynamics lit- 
erature these order-forming transitions are usually referred to as non-equilibrium transitions, to stress the fundamental 
difference in the role of noise with respect to the classical case of equilibrium transitions, which exhibit an increase in 
disorder as the amplitude of internal fluctuations increases. 

Here, we propose an overview of the main stochastic processes related to the presence of Gaussian white noise, 
focusing on the fundamental mechanisms able to induce spatial coherence. We concentrate on Gaussian white noise 
because it provides a reasonable assumption for the unavoidable randomness of real systems - the spatial and temporal 
scales of the Gaussian white noise are much shorter than the characteristic scales over which the spatio-temporal 
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Figure 1: Example of aerial photographs showing vegetation patterns (tiger bush), (a) Somalia (9°20'N, 48°46'E), (b) Niger (13 21'N, 2°5'E), 
(c) Somalia (9°32'N, 49°19'E), (d) Somalia (9°43'N, 49°17'E), (e) Niger (13°24'N, 1°57'E), (f) Somalia (7°41'N, 48°0'E), (g) Senegal (15°6'N, 
15°16'W), and (h) Argentina (54°51'S, 65° 17'W). Google Earth imagery © Google Inc. Used with permission. 
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dynamics of the field variable are evolving - and, therefore, it is typically adopted in stochastic modeling. Moreover, 
the white noise assumption simplifies analytical and numerical calculations. We call "patterned" a field that exhibits an 
ordered state with organized spatial structures. This general definition, including both periodic as well as multiscale 
patterns, is often adopted in the environmental sciences, where the number of different processes can prevent the 
organization of the system with a specific wavelength. We define multiscale those patterns that are scale-free (36] [37), 
in the sense that spatial coherence emerges without showing a clear periodicity. Depending on their behavior in time, 
patterns can be also classified as steady or transient, on the basis of whether the spatial coherence is constant in time 
or appears only temporarily with the tendency to fade out with time. Steady patterns are here defined as statistically 
steady in time. This means that, once the steady state is reached, the field variable can locally assume different values, 
but the mean characteristics of ordered spatial structures remains the same. In the case of oscillating patterns the 
spatial coherence instead fluctuates with time and patterns periodically emerge and disappear. 

In Section [2] a mathematical model of the spatio-temporal dynamics is introduced. We then consider two simple 
stochastic models (Sections[3]and|4]), in order to clarify the interplay among three fundamental mechanisms: the local 
dynamics, the noise component and the spatial coupling. Pattern formation with temporal phase transition is described 
in Section [5] Concluding remarks are given in Section [6] Analytical prognostic tools and numerical algorithms for 



pattern detection are treated in Appendix A 



2. Stochastic Modeling 

The spatio-temporal dynamics of the state variable, (/>, can be expressed, at any point r = (x, y), as the sum of 
four terms: (i) a function, f(<fi), of local dynamics; (ii) a multiplicative noise term, g{(/>)^{r,t); (iii) a term, D £[</>], 
accounting for the spatial interactions with the other points of the domain, and (iv) an additive random component 
f a (r, t). Therefore, the dynamics read 

^ = f(c/>) + g(0#r, f) + D£W + &(r, t), (1) 
at 

where £, is an operator expressing the spatial coupling of the dynamics, while D is the strength of the spatial coupling. 
The description of the spatio-temporal stochastic resonance and coherence - two mechanisms of noise-induced pattern 
formation that need the cooperation of a temporal periodicity - is not included in this review. Thus, we will concentrate 
on the case of dynamics in which the state of the system is determined by one state variable, (p, without time-dependent 
forcing terms. 

The crucial point in the dynamical systems here investigated is that pattern formation is noise-induced, i.e., it is due 
to random fluctuations and does not occur in the deterministic counterpart of the dynamics. In fact, these symmetry- 
breaking states vanish as the noise intensity drops below a critical value depending on the specific spatiotemporal 
stochastic model considered. For some configurations these noise-induced transitions are re-entrant. This means that 
the ordered phase is reached beyond a threshold but is then destroyed if the noise intensity exceeds a higher threshold. 
In these cases, the noise has a constructive effect only when its intensity is within a certain interval of values. Smaller 
or larger values are either too weak or too strong to induce ordered structures. We consider a white (in time and space) 
Gaussian noise with zero mean and correlation given by 

<f (r, f)£(r', f')> = 2*J(r - r')6(t - ?'), (2) 

where s is the noise intensity. We interpret the Langevin equation ([T]i in the Stratonovich sense. In this case 
{g{<p)^{<p)) - s{g{<p)g' '((/>)) , where g'((p) is the derivative of g with respect to <p. In contrast, under Ito interpretation, 
one has <g(0)?(0)> = 0302). 

A number of mathematical models can be used to express the spatial coupling in spatiotemporal dynamics. We 
call pattern-forming those operators that, under suitable conditions, are able to generate periodic patterns even without 
noise. In contrast, non-pattern-forming operators are able to give spatial coherence inducing multiscale patterns, 
without selecting a clear dominant length scale. A typical example of non-pattern-forming operator is the Laplacian, 
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which is widely used to represent the effect of the diffusion mechanisms in a dynamical system. This operator accounts 
for spatial interactions between a point of the domain and its nearest neighbors, and is therefore considered as a short- 
range spatial coupling. 

A mathematical structure able to describe pattern-forming couplings is instead the Swift-Hohenberg operator 

S<H[<f>]=-(V 2 + 1&% (4) 

where Icq is a parameter corresponding to the wavenumber selected by the spatial interactions. It should be noted that, 
beside the effect of short-range interactions expressed by the Laplacian operator V 2 , the biharmonic term, V 4 , accounts 
for long-range interactions. Indeed, in a finite difference discrete representation of equation (|4]) the biharmonic oper- 
ator accounts for interactions with points of the domain located next to the nearest neighbors. The Swift-Hohenberg 
operator is one of the simplest types of coupling able to account for both short and long range interactions and to 
form periodic patterns. For this reason it has been widely adopted in different applications ll24l[32l . The structure Q 
was first introduced by fl] to study the effect of hydrodynamic fluctuations in systems exhibiting Rayleigh-Benard 
convection ||25ll . However, the interplay between short and long range interactions expressed by Q is a recurrent 
mechanism of pattern formation in nature. For example, in landscape ecology, cooperative interactions for vegetation 
growth - such as mulching, shading, absence of biological crusts [38, 39l l40l |4T1 l42ll - occur in the short range of 
plants' crown areas, while inhibitory effects hindering vegetation establishment - such as competition for water and 
nutrients through the root system BT1 [38l l43l l44l - typically occur at larger distances. 



3. Additive noise 

Consider the stochastic model 

^=a0 + ZX£[0]+&(r,t) (5) 
at 

where <p{r, t) is the scalar field, a is a parameter, and £,( r > t) is a zero-mean Gaussian white (in space and time) 
noise with intensity s a . Equation |5]l is the prototype model used to show how patterns may occur in the absence 
of multiplicative noise (i.e., g(<p) — in the general equation ([TJ) and of a time-dependent forcing. We concentrate 
on linear deterministic dynamics to point out the fundamental mechanisms able to induce pattern formation, without 
invoking nonlinearities, which in this case do not substantially change the pattern properties. In so doing, we present 
very simple, though common and realistic, models of pattern formation. In this section we will first study the case 
where £,[(p] is a pattern forming spatial coupling. 

3.1. Pattern fanning coupling 
The prototype model is 

^ = fl 0-D(V 2 +£ 2 ) 2 + £,. (6) 

The deterministic part of the dynamics does not generate patterns for any value of a: if a < the system is damped 
to zero without showing any spatial coherence, if a is positive, no steady states exist and the dynamics of <p diverge 
without displaying any ordered spatial structures. Additive noise, is able to keep the dynamics away from the 
homogenous deterministic steady state even though in the underlying deterministic dynamics /(</>) would tend to 
cause the convergence to the homogenous state. In these conditions patterns emerge and are continuously sustained 
by noise. These patterns are noise-induced in that they disappear and the homogeneous stable state <p — is restored 
if the noise intensity is set to zero. Figure [2] reports some results from numerical simulations, including the spatial 
field, the probability density function (pdf) of <p, and the azimuth-averaged power spectrum and the structure function 
of the stochastic model i5l. More details on the numerical methods are provided in Appendix A The results confirm 



that a very clear and statistically stable pattern occurs in spite of a being negative: the noise component moves the 
dynamics away from the deterministic steady state 0o = and allows the spatial differential terms to drive the field 
into a patterned state with wave length 2n/ko. 
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Figure 2: Model {6} at f = 100, with a = -1, D 
tones are used for positive and negative values of 4 
numerical simulations, dotted: structure function). 



= 10, kg = 1, and s a = 0.5. First panel: numerical simulations of the field. Black and white 
respectively. Second panel: pdf of 0. Third panel: azimuthal-averaged power spectrum (solid: 



This pattern-inducing role of the noise can be detected through the structure function, as defined in Appendix A. 3 



Using Equation ( A. 12 1 one obtains 



D(-k 2 + kl) 2 - : 



(7) 



Thus, even for a < the steady state structure function has a maximum at k = ko (see third row of figure |2j 
where the numerical power spectrum and the structure function are compared at steady state). This result confirms 
that additive random fluctuations are able to induce a stable pattern. When the noise is absent (s a = 0) the steady state 
structure function is uniformly null and no patterns form. 



The normal mode stability analysis (see Appendix A.l i is unable to detect the occurrence of patterns. In fact 



according to this analysis no pattern should emerge when a < 0. Similarly, the generalized mean field technique (see 



Appendix A.4. 1 1 is unable to capture the constructive role of additive noise. 
The classical mean field technique (see | Appendix A. 4~2| ) leads to 



df 



= /W+«-o#- fl 



( 20 



8C 
A 2 



(</>i -m) + tjaj, 



(8) 



where A is the spatial step (see Appendix A I. This analysis shows that the order parameter, m, does not change, i.e., 
the periodic patterns induced by additive noise do not entail phase transitions. The pdfs of the field variable, </>, shown 
in Figure [2] confirm the theoretical findings from the mean field analysis. Indeed, one observes that the pdfs remain 
unimodal and symmetrical at any time, with the mean at <p — in spite of the appearance of patterns. 



3.2. Non-pattern forming coupling 

In this section, we consider the same interplay between the local deterministic component, f {(/>), and the noise 
component investigated in the previous section. However, we consider a spatial coupling which is not able to select a 
specific wavelength. The prototype model becomes 



d<$> 



• + DV 2 (P + t a . 



(9) 



The introduction of the additive random component allows one to obtain very interesting patterned fields. The effect 
of additive noise is even more surprising than with a pattern-forming coupling. In that case, (unsteady) patterns were 
in fact already potentially present in the deterministic dynamics (see Section 3.1 1. In contrast, here the deterministic 
dynamics does not reveal any transient spatial coherence. 

Considering the steady state structure function, one obtains 



S(k,t) = 



Dk 2 



(10) 
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Figure 3: Model J9} with a = —0.1, £* = 2.5, and s a = 5. The columns refers to 0, 200, and 400 time units. First row: numerical simulations of the 
field. Second row: pdfs of (solid: numerical simulation, dotted: classic mean-field, dashed: corrected mean-field). Third row: azimuthal-averaged 
power spectrum (solid: numerical simulations, dotted: structure function). 
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which exhibits a maximum at k — (see third row of figure [3] where the numerical power spectrum and the structure 
function are shown at steady state). No specific periodicity is selected, but a range of wave numbers close to zero 
compete to give rise to multiscale patterns. However, since it is difficult to predict the characteristics of these patterns 
only by looking at the properties of the structure function, we have to rely mostly on numerical simulations. 

An example is reported in figure [3] where results are shown in terms of the spatial field, the pdf, the azimuth- 
averaged power spectrum and the structure function. As expected from the analysis of the structure function, no 
clear periodicity is visible, but many wave lengths are present. The boundaries of the coherence regions are irregular 
and these spatial structures fall then in the class of multiscale fringed patterns, which are especially relevant in the 
environmental sciences [36 , 37 1. In fact, a number of environmental patterns exhibit a spatial behavior very similar to 
the one shown in figure [3] A typical example is the distribution of vegetated sites in semi-arid environments (8). 

In this case, other prognostic tools, such as the modal stability analysis and the generalized mean field theory, fail 
to provide useful indications. The pdf is unimodal (see the second row of the figure [3j and its mean coincides with 
the basic homogeneous stable state (i.e., m — <po = 0). Therefore, there is no phase transition. This behavior follows 
the general rule that Gaussian additive noise is unable to give rise to phase transitions (i.e. changes of m) for any type 
of spatial coupling. 

Pattern formation induced by additive noise is usually introduced in the scientific literature as a remarkable exam- 
ple of noisy precursor near a deterministic pattern-forming bifurcation ll32l . The additive noise acts on a deterministic 
system that exhibits a bifurcation point between a homogeneous stable state and a stable patterned state (an example 
is the Ginzburg-Landau model ll33l [321 ). In this case, the role of the additive noise is to unveil the intrinsic spatial 
periodicity of the deterministic system even before reaching the pattern-forming bifurcation. 

This point of view suggests that a deterministic bifurcation is necessary in order to have an additive noise generating 
a pattern. The example with f{<p) = cup we have just presented demonstrates, instead, that this is not necessarily true. 
In this case there is no bifurcation since the dynamical system diverges when a > 0. Therefore, the existence of a 
deterministic bifurcation is not a necessary condition for pattern formation. Patterns emerge as an effect of additive 
noise, which unveils the capability of the deterministic component of the dynamical system to induce transient periodic 
patterns also when the asymptotic stable state is homogenous. Thus, noise exploits this capability and hampers 
patterns to disappear. 

Moreover, it should be noted that the presence of nonlinear components in the deterministic dynamics does not sub- 
stantially change any of the previous results. The fine details of the patterns can change, but neither their stable 
occurrence nor their dominant wave length (if detectable) changes. 



4. Multiplicative noise 

The cooperation between multiplicative noise and spatial coupling is based on two key actions: (i) the multiplica- 
tive random component temporarily destabilizes the homogeneous stable state, <po, of the underlying deterministic 
dynamics, and (ii) the spatial coupling acts during this instability, thereby generating and stabilizing a pattern. The 
basic model is 

^=/(0 + g<0?(r,f) + aC[0], (11) 
at 

where, with respect to the general equation ([TJ, has been eliminated in order to isolate the role of the multiplicative 
noise. £ is a zero-average Gaussian white noise with intensity s. We indicate with <f> the stable homogeneous state 



of the system in the deterministic case. Namely, <f>(r, t) = (po is a homogeneous solution of ( 1 1 1 when s — (i.e., 
f((po) — because -C[<p] = in homogeneous states). Moreover, we consider cases where g(<f>o) = 0, so that the noise 
does not have the possibility to destabilize the homogeneous steady state. 



The analytical tools detecting the possible presence of the short term instability are described in Appendix A. 2 
For values of s lower than a critical value, s < s c , the state variable <p(x, t) experiences fluctuations about 0o but noise 
does not play any constructive role. The system remains blocked in the disordered phase and no patterns occur. Only 
transiently, the spatial coupling might be able to induce patterns that disappear as the system approaches its steady 
state. Conversely, when the noise increases above a critical level, s > s c , the spatial term can take advantage of the 
noise-induced short term instability and prevent that the displacement from the homogeneous equilibrium state decays 
to zero. In this way, the spatial coupling traps the system in a new ordered state, maintaining the dynamics far from 
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Figure 4: (a) Behavior of (if>) - 0o obtained as ensemble average of 10 6 realizations of the model jl3| with a = — 1, The initial condition is = 0.1, 
s = 0.5 and s = 5 for the dashed and solid curves, respectively, (b) Dispersion relation of the model (12) , with s = 0.1,1,2 (dotted, dashed and 
solid curves, respectively), a = -1, D — 1, and Icq = 1. 



the state <pQ. Equation ( 1 1 1 is here interpreted in the Stratonovich sense, where (g(<p)^((p)} = s(g((p)g'(ip)}, while under 



Ito's interpretation no short-term instability occurs, as (g((p)^(<p)) = (see Appendix A.2\ 



Other works have been proposed to describe the role of multiplicative noise in pattern formation and phase tran- 
sitions Il34l [35l l45l . and to show how multiplicative noise can induce periodic patterns ll46l l33l . However, they all 
present quite complicate non linear expressions to represent the local dynamics and the noise terms, so that their 
physical interpretation is not always straightforward. For example, when g(<po) + 0, the noise plays a role similar to 
£,(r, t) in Eq. ([5]). Moreover, the results of those models are qualitatively similar to those described in the following 
sections, provided that the interplay between short-term instability and spatial coupling remains the same. 

4.1. Pattern forming coupling 

To illustrate how pattern formation can be driven by multiplicative noise in the presence of a pattern-forming 
spacial coupling, we consider the model 

^ =a0-0 3 +0<f-D(£5 + V 2 ) 2 0, (12) 

where a is a negative number, the random component is modulated by a function g((p) = <p, and £ is a zero mean 
white Gaussian noise, with intensity s. The local dynamics are f(<p) = a<p - 3 , where the nonlinear term -0 3 has 
been introduced to avoid that the dynamics diverge. Indeed, the linear term prevents the dynamics from relaxing 
to a homogeneous steady state, thereby allowing the spatial terms to be different from zero, while the nonlinear 
term prevents the divergence of the dynamics (i.e., it ensures the convergence to a statistically stable state). The 
deterministic homogeneous stable state, obtained as a solution of f((po) = 0, is (f>{r, t) = O = 0. From the short-term 



instability (see Appendix A. 2 1 one has that <po is stable for s < -a and becomes instable for s > -a. To show this 
point, Figure reports the time behavior of the ensemble average of a number of numerically evaluated realizations 



of the zero-dimensional stochastic model obtained eliminating the spatial component from equation ( 12 1, that is 

^7 = f{<P) + 8(mt) = a<f>-</, 3 + tf(t). (13) 
at 

It is evident from figure [4^ that the growth phase appears only when s > s c = -a and at the beginning of simula- 
tions (i.e., at short term), while the effect of the initial perturbation disappears in the long run. Mathematically, the 
short-term instability can be understood by observing that when <p is close to zero the disturbance effect due to the 
(multiplicative) noise tends to prevail on the restoring effect of /. When <p grows, the leading term, tp 3 , prevails on g 
and the local dynamics, /, tend to restore the state 0r> 
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Figure 5: Model (12) aU = 100, with a = -1, D = 15, /t = 1, 
Third panel: azimuthal-averaged power spectrum. 



2.5. First panel: numerical simulations of the field. Second panel: pdf of $ 



Once the presence of a short-term instability has been detected, the capability of the spatiotemporal stochastic 



model ( 12 1 to give rise to patterns can be investigated through the stability analysis by normal modes, see Appendix 



A. 1 The dispersion relation is 



y(k) = a + s- D(Icq 



k 2 ) 2 , 



(14) 



which provides the same threshold s c for the neutral stability found with the short-term instability, while the maximum 
amplification is for the wavenumber k = ko (see figure^). It follows that statistically steady periodic patterns, with 
wave length A = 2n/ko, emerge when the noise intensity exceeds the threshold s c = -a. Moreover, the critical value 



of the noise intensity for the neutral stability, s c = [-a + Z)(&o - k ) ], is confirmed by the structure function (see 
|Appendix A.3| >. The generalized mean-field analysis of the most unstable mode predicts an unconfined region of 
instability. This means that, if s > s c , the system is able to exhibit periodic patterns for any value of D. 

The numerical simulation of the stochastic model ( 12 1 confirm these theoretical findings. Figure [5] shows an ex- 
ample of patterns emerging from the simulation. These patterns have the same basic characteristics as those observed 
in the case of additive noise (see figure |2j. They are statistically stable and exhibit a clear dominant wave length cor- 
responding to ko- In this case the pdf of the field is weakly bimodal with the zero mean, demonstrating that no phase 
transition occurs. However, there are some important differences with respect to the case of additive noise. Firstly, the 
boundaries appear to be more regular when patterns are induced by multiplicative noise. Such aspect is also displayed 
by the power spectrum of the field, which shows a more sharp peak at ko in the case of multiplicative noise (compare 
the third rows of figures[2]and[5]). This difference is due to the fact that multiplicative noise is modulated by the local 
value of <p, and this has the effect to make the boundaries of the pattern more regular because the (p field is spatially 
correlated (from the definition itself of patterned state). Another difference is that patterns induced by additive noise 
exhibit more stable shapes than those emerging as an effect of multiplicative noise. For example, patterns shown 
in figure [5] seem to evolve from a labyrinthine shape to a striped shape. Overall, numerical simulations show that 
patterns induced by additive noise quickly reach their steady state without showing any transient temporal evolution, 
while those induced by multiplicative noise present a transient behavior, during which they modify their shape until 
the steady state is reached. The third difference is the possible occurrence of a weak bimodality in the pdf of <p in the 
dynamics driven by multiplicative noise (see figure[5]). The presence of such bimodality depends on model structure, 
parameter values, and field size; however, it generally remains weak. 

Let's now look at the effect of nonlinear g(<p) terms on the short-term instability. Consider equation ( 12 1 with 
a < and g{(j)) = <p a . If a > 1 no short-term instability occurs. This result is explained interpreting the short-term 
behavior as a balance between the tendency of f(<f>) to restore the homogeneous state, (f> — (/>q = 0, and the diverging 
action of g(<f>) ■ Since (f> is close to zero, the power a > 1 of the function g{<p) reduces the effect of the noise term, 



which becomes unable to contrast the action of the leading term, a<p, of f((f>). Patterns occur only transiently and the 
field then rapidly decays to the homogeneous state cpo. 

Conversely, if a < 1, short-term instability occurs in that the multiplicative component always overcomes the action 
of f((p) close to zero. The balance between f(<f>) and g((p)^ is reverted when <p moves away from 0o = 0, and this fact 
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hampers the dynamical system to diverge. In this case, patterns are statistically stable and exhibit the same dominant 
wave length as those shown in figure [5] 



4.2. Non-pattern forming coupling 

We explore the capability of spatio-temporal models driven by multiplicative noise to generate patterns when a 
non-pattern forming spatial coupling is adopted. To focus on the role of the type of spatial operator, we consider the 



same model as in Section 4. 1 but with a diffusive Laplacian operator, namely 

^ =CKt>-4> i + 4>% + DV 2 <p. (15) 
where £ is a zero mean white Gaussian noise, with intensity s, and equation ( fl5] > is interpreted in the Stratonovich 



sense. The model ( 15 1 has a homogeneous deterministic stable state at <po = and exhibits short-term instability when 



s > s c = -a (see Section 4.1 1. The dispersion relation obtained with the linear stability analysis reads 



y{k) = a + s-Dk 2 . (16) 

We can make the following three remarks. First, the value s = -a of the noise intensity marks the condition of marginal 
stability. No unstable wave numbers occur when s < -a, while the wave numbers lower than V(,s + a)/D become 
unstable if s > -a. The threshold s = -a coincides with the one obtained in the short term analysis. Second, the 
strength, D, of the spatial diffusive coupling impacts the range of unstable wavenumbers. In particular, the unstable 
wave numbers decrease when D increases, consistently with the fact that the diffusive coupling introduces spatial 
coherence in the random field. However, D does not impact the occurrence of instability, in that it depends only on 
the noise intensity. Third, the most (linearly) unstable mode is always k max = 0, for any D and provided that s > -a. 
The classic mean-field analysis predicts a phase transition for any value of the spatial coupling, D. In other words, 
once the critical threshold of noise intensity is exceeded, the system is always able to move to a new ordered state. 

Figure|6]shows that the emergence of a pattern with no clear periodicity. It evolves in time and tends to disappear 
in the long term. The pdf of the field reveals that patterns occur during a phase transition from the initial basic state 
having order parameter m = 4>o — to a new substantially homogeneous state with m + 0. In particular, for the case 
shown in figure |6]numerical simulations give m = 0.8 for t > 100 time units. Although with diffusive coupling both 
the cases of additive and multiplicative noise (not shown for sake of brevity) exhibit a well defined peak at k — 0, 
numerical simulations show a different scenario (compare figures [6] and [3j. While the presence of additive noise 
produces steady multiscale fringed patterns, multiplicative noise induces only smooth transient patterns. The reason 
of this different temporal behavior is that, in the multiplicative case, the diffusive operator is unable to maintain the 
system far from homogeneous condition, in spite of the initial instability. Indeed, to have steady patterns sustained 
by multiplicative noise the presence of a pattern-forming spatial coupling is necessary. When other types of spatial 
couplings are considered, they are unable to block the system far from the homogeneous state. In these cases, the 
spatial coupling interacts with the short-term instability, as detected by the dispersion relation - recall that the stability 
analysis is performed on an equation that approximates only the first stages of the ensemble average dynamics - but 
this interaction lasts only until the temporal dynamics are able to sustain the instability. Thereafter, the patterns 
undergo the same fate as the initial instability, i.e., they tend to disappear. In the long run the main legacy of the 
spatial coupling is the phase transition (i.e., m + 0), though for a homogeneous field. Similarly to the case discussed 



in Section 4. 1 the coherence regions are much smoother than those observed with additive noise. Even in the this 
case, this difference is due to the multiplicative nature of the noise, which entails that the g((p)^ term is spatially 
correlated. 



5. Patterns with temporal phase transition 

In this section we consider the extension to spatial systems of noise-induced transitions in purely temporal systems 
BTll . Such transitions correspond to the occurrence of steady state probability distributions whose modes are different 
from the equilibrium states of the corresponding deterministic system. A relevant case is represented by systems 
exhibiting noise-induced bistability. In this case, suitable noise intensities are able to generate pdfs with two modes 
even though the deterministic dynamics have only one stable state. Two key ingredients are needed to activate this 
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p(<A) p(0 P(0) 




Figure 6: Model <JT3J with a = — 1, s = 2, D = 5. The columns refer to 0, 10, and 40 time units. First row: numerical simulations of the field. 
Second row: pdfs of 0. 

type of stochastic dynamics. First, a deterministic local kinetics, f(<p), which tends to drive the dynamical system 
towards the steady state, (f> — <pQ. Second, a multiplicative random component that tends to drive the state of the 
system away from cp = fa; the intensity of this component is generally maximum at <p — <pQ. As a result of the balance 
between deterministic and stochastic components, bimodal probability distributions of (f> may emerge at steady state. 

In spatiotemporal dynamical systems the spatial coupling could (i) cooperate with the stochastic component to 
prevent the relaxation imposed by the local dynamics and maintaining the system away from the uniform state, {(p) = 
<po, and (ii) give spatial coherence to the field creating a patterned state where the coherent regions correspond to the 
two modes existing in the underlying temporal dynamics. 

The mechanism here discussed is sometimes called entropy-driven pattern formation [32 1 as the dynamical system 
escapes from the minimum of the potential (i.e., (f> = <po) because of the strength of noise (which is an entropy source). 
There are two main differences with respect to the case presented in the previous sections: (i) patterns do not result 
from a short-term instability, and (ii) they emerge even if the noise is interpreted according to Ito's rule. 

5.1. Model with g(<p ) = 

It is interesting to consider what happens when no noise term is present for <p = <pQ, i.e. when g{<po) = 0. Indeed, in 
this case the noise component is unable to unlock the system from the deterministic stable state <po and to sustain the 
pattern-forming effect of the spatial coupling by maintaining the dynamics away from the homogeneous stable state, 
<p = <p . 

We consider the model 

^ = -a0 + <p(l - mt) - D(kl + V 2 f<p, (17) 

interpreted according to Ito. As in the case of section [3] we concentrate on the case of linear local dynamics to show 
how patterns may emerge even without invoking nonlinearities in the underlying deterministic dynamics. The purely 
temporal version of the model ( [17} (i.e., d(/>/dt = —a<p + <p(l - <p)t;) shows a noise-induced transition for s c = 4a. The 
steady state pdf reads 
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Figure 7: Model {T7} under Ito interpretation at t = 100. The parameters are a = 0.001, D = 10, and s = 1. 



P(<P) 



(1 



\--2 



1+2 



-Exp 



s(l 



•6]0,1[ 



(18) 



and has always a mode for <p — > 0, while a second mode occurs when s > s c . The onset of bimodality in the model ( fTT) 
is due to cooperation between the noise and the natural boundaries at — and = 1. When the noise is sufficiently 
strong, the system tends to move away from = 0, but the boundary at = 1 prevents the system from visiting the 
whole real axis, and an accumulation of probability close to the upper limit of the domain emerges. 



Since the noise is interpreted according to Ito, model ( 17 1 does not present any noise-induced short term instability 
(i.e., {g(<f>)<;} = 0). However, the spatial coupling is able to exploit the temporal noise-induced transition and to show 
patterns. In figure[7]an example is reported. In spite of the pdf of the temporal system displaying a strong bimodality 
for high enough noise strength, the pdf of the field is unimodal and centered in zero. The generalized mean-field 
analysis is not able to provide further information. 



Patterns can emerge also when in equation ( 17 i the diffusive spatial coupling is used in place of the Swift-Hohenberg 
operator, 



c90 



-a0 + 0(l-0)^(f) 



• DV 2 0- 



(19) 



Figure [8] shows an example of these patterns. For high enough values of the ratio D/a, the pdf is unimodal and 
converges to = 1 . However, the classic mean-field analysis is not able to capture any phase transition of the system. 



6. Conclusions 

We presented different stochastic mechanisms of spatial pattern formation. They all describe spatial coherence 
and organization as noise-induced phenomena, in the sense that these patterns emerge as an effect of the randomness 
of the systems drivers. 

Additive noise plays a fundamental role when the deterministic local dynamics tends to drive the field variable 
towards a uniform steady state, while noise is able to maintain the dynamics away from the uniform steady state. The 
interaction of additive noise with the spatial coupling provides a simple and realistic, mechanism of pattern formation. 
In the presence of a multiplicative component of adequate intensity the spatial coupling exploits the initial instability 
of the system to generate ordered structures, which in the absence of noise would tend to disappear in the long run. 
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Figure 8: Model j!9[ at t = 300 under Ito interpretation. The initial conditions are given by uniformly distributed random numbers between [0.49, 
0.51]. Black and white tones are used for the value intervals [0.5, 1] and [0, 0.5], respectively. The parameters are a = 0.001, D = 25, s = 40. 



The stochastic models presented here show how noise may play a crucial role in pattern formation. However, 
most of the literature on self-organized morphogenesis in the environment is based on deterministic mechanisms. The 
limited application of stochastic theories to environmental patterns is likely due to the fact that most of the stochastic 
models use some specific (and complicated) non linear terms both in the local deterministic dynamics and in the 
multiplicative function, g((f>), of the noise component. The use of these "ad hoc" functions limits the applicability 
of these theories to process based environmental modeling. Thus, only few studies have investigated the possible 
emergence of vegetation patterns as a noise-induced effect. However, because noisy fluctuations - such as those 
associated with fires, rain, or soil heterogeneity - are a recurrent feature of environmental drivers, their randomness 
can actually induce spatial coherence in a number of environmental processes, as well as in problems related to front 
propagation under non-equilibrium conditions 1 48 , 49 32] . 



Appendix A. Analytical and numerical tools 

The mathematical complexity of the spatio-temporal models of type ([T]i hampers general analytical solutions. For 
this reason, several approximate analytical techniques have been developed in order to obtain some indications of 
pattern formation. We briefly describe the most important ones. 

Appendix A.l. Linear stability analysis by normal modes 

When the occurrence of a dominant wavelength is the main symptom of pattern formation, the first available 
prognostic tool is provided by the normal mode linear stability analysis. This analysis is based on the idea of disturbing 
the basic state of the system with a hypothetical infinitesimal perturbation, and to assess whether the perturbation 
grows in time (in which case patterns have the possibility to emerge) or not. The analysis involves three steps. 
Firstly, a deterministic equation for the spatiotemporal dynamics of the ensemble average of the field variable, (<p), 
is determined and its homogeneous steady state is found. In general, the multiplicative random component can be 
expressed as (g(0)£) = s(gs(<f>)), where gs((f>) is a function of the state variable. Applying Novikov's theorem [33 1, 
for the case of Gaussian white noise interpreted in the Stratonovich sense, we have (g(4>)^) = s(g((p)g'(<p)), i.e. 
gs(<P) = gWg'W- Instead, with Ito's interpretation we have (g((/>)0 - 0. Using equation ([TJ, we find 

= </(*)> + s{g s {<l>)) + D£[(cf>}]. (A.l) 

at 

The basic state, {(f>) = </>q, is obtained as the zero of equation ( |A. 1) > at steady state, i.e. f(4>o) + s gs(</>o) - 0- 
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Equation (A.l i is then linearized, and Taylor's expansion of equation (A.l I around <p = <Po, truncated to the first 
order, provides 



d{<t>) 
dt 



f'((l>o)(<p) + sg' s (cf> )(<l>)+D£[(cf>)], 



(A.2) 



where f'(</> ) = 



dJU/A 



and g' s (cf>o) = 



The basic state (0) = <po is perturbed (third step) by adding an 



infinitesimal harmonic perturbation 



(cf>) = 0o + 0e y,+/kr , (A.3) 
with (j) being the perturbation amplitude, y the growth factor, i = V— 1 the imaginary unit, k = (k x , ky) the wave 



number vector of the perturbation, and r = (x,y) the coordinate vector. If equation (A.3i is inserted in (A.2i, one 
obtains the so-called dispersion relation 

y(k) = /'(0 O ) + sg'sifo) + Dh L (k), (A.4) 
where li£(k) is a function of the wave number k = |k| which depends on the specific form of spatial coupling X. 



considered. From equation (A.4i the threshold value of the noise intensity is easily obtained by setting the marginal 
condition, y — 



f'(<p a ) + Dh £ (k) 
Sc g' s (4>o) ' 

c When s > s c , the growth factor, y, is positive and spatial patterns may occur. 



(A.5) 



Appendix A.2. Short-term instability 

A second tool can be used to assess the possible short term instability in the dynamical system. The transient 
instability is important because it tends to move the dynamics away from the basic state (<p) - <pQ. If this phenomenon 
is accompanied by a suitable spatial coupling, the system can be trapped in a new ordered state. The first steps of the 



stability analysis are the same as those described in Appendix A. 1 and lead to Eq. ( A. 1 1. Since we are here interested 
in the initial evolution of small displacements from <f>o, we can assume (/(</>)) ~ f((<p)) and (gs(4>)) ~ gs ((</>))■ For 
the same reason, we can also neglect the spatial gradients of the fluctuations assuming that in the short term they are 
small. Equation dATh can be approximated at any point of the field as ll32l 



dt 



/«0» + iW«0» = /«//«0». 



(A.6) 



where f e jy is often indicated as the effective kinetics. The short term stability analysis of the state (po by equation 
(A.6 1 concerns only the temporal dynamics at a generic point of the field. 

Equation (A.6 1 clearly shows that noise can destabilize the deterministically stable state <pQ. This can happen in 
two possible ways: (i) when the roots of f e ff((<p)) = do not coincide with those of /((</>)) — or (ii) when the state 
<po remains a zero of the r.h.s. of equation ( |A.6[ ), but sufficiently high noise intensities destabilize this state an unstable 
one; this occurs when the following condition is met 



d feff 



d(4>) 



> 0. 



(A.7) 



In both cases noise has to be multiplicative for a transition to occur. The noise threshold s c can be obtained by setting 



the inequality in (A.7 1 equal to zero. The condition (|A.7|i derives from the first-order truncated Taylor expansion of 



the function f e ft around <f> , which yields 



d(<P) ^ df e ff 

dt ~ d{(f>) 



((/>)■ 



(A.8) 
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When f e ff((f>o) - 0, the sign of the coefficient of (<p) on the r.h.s. of the previous relation determines the stability of 
small perturbations around <p^. 

When the noise is additive the effective kinetics are f e ff((<p)) = /(((/>)). Thus, the stable states of (A.61 are the 
same as those of the deterministic counterpart of the process. 

It is also worth to stressing the impact of the type of noise and of its interpretation. While in the Stratonovich case, 
gs(4>) = sg(<f>)g'(<l>), using Ito interpretation we have gs((f>) = 0. Therefore, the noise-induced instability is possible 
only in the Stratonovich interpretation of the Langevin equation ([T]) and it cannot occur when Ito's framework is 
adopted. 



Appendix A. 3. Structure function 

The presence of patterns modifies the correlation structure of the field. Instead of considering the correlation 
function, this method analyzes its Fourier transform in space, which is known with the name of structure function and 
defined as S(k, t) = {<p(k, t)$>(-k, f)), where k, f) is the Fourier transform of <p(r,t) and k = (k x ,k y ) is the wave 
number vector. The structure function is therefore equal to the power spectrum of the field <p. The first-order temporal 
derivative of the structure function reads 



8S(k,t) Id4>(k,t) 



at 



dt 



-0(-k,O) + 



d4>(-k, t) 
dt 



<p(k,t) 



(A.9) 



By taking the Fourier transform of the linearized version of Eq. ([T]i (with F(t) = 0), which is obtained by following 
similar steps to those leading to equation (A. 2 1, and considering the noise terms as white in space, we obtain 



<50(k, f) 
dt 



f'iMkk, t) + g'(<fa)fck, m) + Ut) + Dh £ (k)fck, t), 



(A. 10) 



where li£(k) is the same operator already defined in Appendix A.l By substituting Eq. ( A. 10 1 into Eq. (A. 9 1 and 
using Novikov's theorem ll33ll to express the terms (^>(k, f)0(-k, f)£( f )) = sS(k, t) and {^>{±k, f)fa(f)) = s a (where s a 
is the intensity of the additive white Gaussian noise), we obtain 



8S(k,t) 
dt 



2 [/* (0„) + Dh £ (k) + g'{cf> )s]S(k, t) + 2s a 



At steady state the structure function reads 

S st (k) = - 



[f'(M + Dh £ (k)+g'(<f>o)s] 



(A. 11) 



(A. 12) 



Eq. (A. 12 1 can be investigated to understand if periodic patterns, corresponding to a maximum of the structure 
function for wave numbers k different from zero, are expected to appear. Equation ( |A. 12) > shows that the additive 
noise is fundamental to have a non-null steady-state structure function. 



Appendix A.4. Mean-field analysis 

The mean-field theory is typically used to provide an approximated solution of stochastic partial differential equa- 
tions for spatially-extended systems. The method is valuable mainly for a qualitative analysis of (stochastic) spa- 
tiotemporal dynamics B4l [351 l50l ISTIl . 

The mean field technique adopts a finite difference representation of the stochastic spatiotemporal dynamics ([TJ 

^ = fm + g(<Pd ■ Ut) + D ■ l{<Pu ft) + h(<Pi) ■ F(t) + UM, (A. 13) 

where </>,-, and , are the values of <p, and at site 2, respectively, i runs across all the cells of the discretized 
domain, and j e nn(i) refers to the neighbors of the z'-th site involved in the discretized representation, Z(0, , <f>j), of the 
specific spatial coupling considered. A general expression for /(</>,, <pj) is 



l(4>i,<f>j) = w4i + 2^ w j<f>j> 

jeim(i) 
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(A. 14) 



where the number of neighbors, nn(i), and the weighting factors w, and Wj depend on the specific finite difference 
scheme adopted to numerically approximate the spatial operator (52). The analytical solution of equation ( A. 1 3| > is 
hampered by the fact that the dynamics of <pj are coupled to that of the neighboring points. To circumvent this issue, 
the mean field approach assumes that (i) the variables <pj can be approximated by their local ensemble mean, (</>j), and 
(ii) there is a relation linking {(pj) to the ensemble average, ((p). 



Appendix A.4.1. Generalized mean-field theory 

To study the stability of the homogeneous steady state with respect to periodic patterns, the pattern is approximated 
by a harmonic function, 



<<^> = <0>COS[k -(I"; -I";)], 



where k = (k x , k y ) is the wave number vector. The function l((f>j, <f>j) in equation (A. 13 1 is approximated as 



l(<pi,(pj) * l h (<f>i,(4>),k x ,ky) 
where 40) depends on the spatial coupling considered. 



(A. 15) 



(A. 16) 



Under the assumption (A. 16 1, the dynamics (A. 13 1 of 0, become independent of those of the neighboring points. 



Thus, it is possible to determine exact expressions for the steady-state probability distributions p s ,(4>; (4>), k x , k y ) of ( 
The self-consistency condition 



(<P) 



oo 



ky) 



(A. 17) 



can be used to obtain the unknown (<p) as a function of k x and ky. 

The occurrence of solutions of equation ( |A.17[ ) different from {(f>) = <po (where <pQ is the uniform steady state of 
the system) corresponds to the loss of stability of the uniform steady state with respect to periodic perturbations. It is 
expected that this loss of stability takes place only for some specific value of the wave numbers k x and ky. 

Other noise-induced phenomena can be investigated with this method. Indeed, spatial pattern formation is not the 
only interesting noise-induced effect. In fact, modifications of other statistical descriptors of the field can be relevant, 
too. Modifications of order parameters are known as phase transitions. In particular, when the spatio-temporal 
average, m, of the state variable at steady state is different from the homogeneous steady state value, a phase transition 
occurs. Non-equilibrium phase transitions are induced by the random forcing. The occurrence of non-equilibrium 
phase transition is neither a necessary nor a sufficient condition for noise-induced pattern formation. Non-equilibrium 
phase transitions imply that noise is able to change the value of the order parameter, but not that ordered geometrical 
structures necessarily emerge. Conversely, we have shown that noise-induced patterns may emerge even when m 
remains unchanged with respect to the disordered case (i.e., no phase transition occur). 



Appendix A.4.2. Classic mean-field theory 

The classic mean field theory can be presented as a simplified version of the generalized mean field. More specif- 
ically, it is assumed that all cells have the same mean, which coincides with the spatiotemporal mean of the field. 



Namely, k = in equation (A. 15 I, i.e. (<pj) = (<p) = m. In this case, equation (A.17i becomes 



<0> 



r*+oo 

4>p st ((/>;(<f>))d4> = F((<f>j). 

<J — oo 



(A. 18) 



The change in the number of solutions of equation ( A. 1 8 1 indicates the existence of a phase transition. The focus of this 
analysis is not on the appearance of periodic patterns but only on the occurrence of phase transitions. The effectiveness 
of this standard mean field approximation can be improved by expressing the values of <f>j in the neighborhood of point 
i as the average between the spatiotemporal mean and the local value of <f> at point i, namely <pj « l/2((<p) + <pj). This 
correction of the mean field approximation accounts for the dependence of <pj on the local conditions |32|. 

The analytical tools here described provide in general some insights into pattern formation, but the definitive way 
to study the noise-induced pattern formation is to numerically simulate the dynamics and assess the emergence of 
spatial coherence through a comparison with homogeneous or disordered states of the system. In fact, the linear 
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stability analysis by normal modes and the short-term instability analysis tend to fail when the noise is additive (i.e. 
g — in Eq. (|T]i). In this case, both analyses predict stable configurations for any noise intensity and any strength of 
the spatial coupling. On the other hand, when only multiplicative noise is present (i.e. £ a = in Eq. ([T}), the structure 
function does not provide any information on pattern formation, as (f> tends to remain equal to zero. 

The typical numerical approach is to discretize the continuous spatial domain using a regular Cartesian lattice 
with spacing Ax = Ay = A. Here we consider a two-dimensional square lattice with 128x128 sites and A = 1. The 
original stochastic partial differential equation ([T]i is then transformed into a system of coupled stochastic ordinary 
differential equations as in Eq. ( [A.13) , Infinitely vast random fields are generally approximated numerically in a 
satisfactory way by periodic boundary conditions, which have been used also in this study. Moreover, unless it was 
otherwise specified, the initial conditions used in the simulations were uniformly distributed random numbers between 
[-0.01, 0.01]. Numerical simulations were carried out with the Heun's predictor-corrector scheme [35],[32|. The pdf 
is obtained from the spatial distribution of <p at fixed time and is numerically evaluated at 100 equally spaced intervals, 
A<f>, that cover the range of (f> values. 
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